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Abstract. The dynamics of homogeneous Robertson- Walker cosmological mod- 
els with a self-interacting scalar field source is examined here in full generality, re- 
quiring only the scalar field potential to be bounded from below and divergent when 



' the field diverges. In this way we are able to give a unified treatment of all the 

already studied cases - such as positive potentials which exhibit asymptotically poly- 
O |' nomial or exponential behaviors - together with its extension to a much wider set of 

physically sensible potentials. Since the set includes potentials with negative inferior 
bound, we are able to give, in particular, the analysis of the asymptotically anti De 
Sitter states for such cosmologies. 



1. Introduction 



Scalar fields have attracted a great deal of attention in cosmology, since the dis- 
covery that such a field can act as an "effective" cosmological constant in driving 



t^. ■ inflation Q18Q. From this first break-trough, which aroused from the simplest con- 
ceivable model, i.e. that of a non self-interacting field, the research on scalar field 
spacetimes has been constantly growing, and quite more general scenarios have been 



considered, such as scalar fields coupled with perfect fluids, non-minimal couplings, 



O ■ and alternate theories of gravity [JTJ, 0, gj, g, 0, g] 

A scalar field spacetime can be viewed, like other matter models coupled with grav- 
ity such as perfect fluids, as a solution of the Einstein field equations which depends 
on the choice of the equation of state for the matter. In the case of the scalar field, 
^ ■ the role of the equation of state is played by the self-interaction potential V , which is 
equal to zero in the "standard" inflationary solution. For several reasons, however, it is 
difficult to believe that this function vanishes. For instance, dimensional reduction of 
fundamental theories to four dimensions typically gives rise to self-interacting scalar 
fields with exponential potentials, coupled to four-dimensional gravity. Therefore ex- 
ponential potentials have been widely considered in the literature (see [ETH and refer- 



ences therein) also with the aim of uncovering possible large-scale observable effects 
(nnL |Z5| , |Zo] , Of course, however, at the present status of our knowledge the spe- 



cific functional form of V is quite unsure and, as a consequence, it would be optimal 
to classify the dynamical behavior of the spacetimes in dependance of all the possi- 
ble choices of the potential function. Relevant attempts have been made also in this 
more general direction. Actually, the application of dynamical systems techniques has 
proven very useful and, due mainly to the works [ lOj |19| , [20| , [21] , |22| , [23| , j24[ ], we know 



the dynamical behavior of scalar field spacetimes for a wide class of non-negative 
potentials; for such potentials therefore, as will be explained below, the contribution 
of the present paper relies in a simplification and completion of known results. 

Non negativity of the potential means that the energy density of the scalar field has 
a positive lower bound. However, in many issues and especially when string theory 
comes into play, it becomes relevant to inspect spacetimes in which the potential still 

l 
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has a lower bound but this bound is negative. From the physical point of view, it 
should be noted that, although care must be given to the fact that local positivity of 
energy density might be violated in these spacetimes, under certain somewhat general 
conditions potentials of this kind generate solutions with positive total energy JT7|]. 
In particular, of course, a constant negative potential generates an "equilibrium" state 
which is just the Anti de Sitter solution (AdS); therefore, a potential with a negative 
minimum generates a class of spacetimes which have an AdS "equilibrium" (see e.g. 
[0. [j~5| [RJ). In the present paper we thus study the dynamical behavior of FRW scalar 
field cosmological models, imposing only very general conditions on the scalar field 
potential; our hypotheses essentially reduce to ask the potential to be bounded from 
below and divergent when the field diverges. Within this quite general framework, 
we give a unified treatment of all the relevant known cases - such as asymptotically 
polynomial and exponential behaviors - as well as its completion to a much wider 
set of possible, physically sensible potentials, also in presence of negative (i.e. AdS) 
minima. It is worth mentioning that our treatment can be applied to the "reverse" 
problem as well, i.e. to homogeneous scalar field collapse. This issue is treated in 
companion papers [O, O, M . 



The paper is organized as follows. First, we formulate the field equations as a 
first order, regular dynamical system for the scalar field, its time derivative, and the 
scale factor. Then, we study in full generality the behavior of the trajectories of this 
system in dependance of the properties of the potential function. Finally, the issue 
of recollapse and that of stability of AdS spacetimes are addressed. In particular, 
we show that that anti deSitter solutions are the only solutions such that <p goes to a 
critical point of V{4>) with negative critical value. 

2. FORMULATION OF THE PROBLEM 

Einstein's field equations for the Robertson-Walker model 

dr 2 + r 2 (d6 2 + sin 2 6 dcp 2 ) 



a) ds 2 = -dt* + a( t y 

coupled with the stress-energy tensor 

(2) = d^d u <P - [^g^d^d^ + g^, 
are given by 

(3) (G° = 8vrT °): -^±^1 = -(0 2 + 2V(0)), 

(X 

(4) (G\ = 8nTl): - {k + 2 ™ = (0 2 - 2V{cj>)). 
in the unknown functions a(t), 4>{t). These equations imply 



(5) = -2<p U + V'(tf>) + 3-<pj = 0. 

We transform the above second order system ©-(|]) into a first order system intro- 
ducing the additional variables h — -,v — We must exclude solutions such that 
4>{t) = 0o constant on some interval but V'(0o) 7^ 0(3> so that we are sure that @ 



khese solutions are easily found by integration of that becomes 
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holds due to the vanishing of the term in brackets. Therefore ((f>,v,h) solves the 
regular system 

(6) 4> = v, 

(7) v = -V'{<f))-3hv 

(8) h = -h 2 -l(v 2 -V^)). 

To see that the converse also holds true, we define the function 

(9) W( ( p,v,h) = h 2 -^(v 2 + 2V(4>)), 

and argue as follows. Given a solution ((f>(t),v(t),h(t)) of @-(|8|), then W(t), the 
derivative of W(t) along the solution, satisfies the differential equation 

(10) W = -2hW, 

that yields W{t) = W exp(— 2 f* h(r)dr). Moreover, we set 

(11) a(t) = a e^ t/l(r)dT , 

where a = |W | -1/2 if W ^ 0, and a(0) = 1 if W = 0. Then the identity 

(12) w ® = ~Wf 

holds, where k = — sgn(Wo), and so <f>(t) and a(t) solve from which @) 

follows. 

We stress some features of the above approach: first, the system obtained is regu- 
lar, with equilibrium points given by (0*, 0, h%) such that V'{4>*) = 0, h 2 = 
(notice that W{<j>*., 0, hj) = 0), that physically correspond to deSitter universe. Given 
a solution ((j)(t),v(t), h(t)), then (</>(—£), —v(—t), —h(—t)) is the time reversed so- 
lution. Moreover, the system represents solutions of all Robertson-Walker cosmolo- 
gies: indeed, from (Q), the set W~ x (0) is invariant by the flow, and by local uniqueness 
this ensures that the sign of W is invariant along the flow; from ( |T2| ) we deduce that 
solutions with W positive, null, or negative, represent scalar field cosmologies with 
k = —1, 0, 1 respectively. From the sign of h(t) one can tell whether, at the instant t, 
the solution is collapsing (h < 0) or expanding (h > 0). 

3. QUALITATIVE ANALYSIS OF THE EXPANDING MODELS 

In this section we consider initially expanding spacetimes, i.e. solutions such that 
h(0) = ho > 0. We suppose that the potential V(4>) is a C 2 function, such that 
limui^oo V(4>) = +oo. We will suppose that critical points of V are isolated, and 
they are either minimum points or nondegenerate maximum points. We also impose 
the weak energy condition to be satisfied at initial time, i.e. Vq + 2V((po) > 0. Of 
course, this is not sufficient for the w.e.c. to be satisfied along the evolution, since 
V((j>) can attain also negative values. It is useful to introduce the energy function 

e = v 2 + 2V((f)), 

in such a way that (§)-(§) imply 

(13) e = -3hv 2 . 

If we let the solution evolve, either there exists T > such that h(T) = 0, or h(t) > 
for all t < sup I, where I is the maximal interval of definition of the solution. In the 
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FIGURE 1. Phase portrait sketch when V(4>) = 5 — 6a; 2 + a; 3 + 2x 4 . 
The behavior when t > for some curves with initial data in the open 
half-space {h > 0} is represented. On the left, the situation nearby a 
local minimum 0* for V, with positive critical value. The three curves 
on the left side, living in the regions where W is positive, negative 
and null, respectively, all converge to the equilibrium point (0*, 0, h*), 
and therefore the corresponding cosmologies are expanding forever. 
On the right, the situation nearby a local minimum of V with negative 
critical value. The two curves, such that W is positive and null, respec- 
tively, leave the half space {h > 0}, therefore passing from expansion 
to collapsing phase. 

latter case, from ( |K)| ) we deduce that |W(i)| < \Wo\ and moreover, recalling (|T3|), 
e < eo = v o + 2V((fio)- Notice also that e(t) is bounded from below since V(0) is, 
and this implies that both <p{t) and v(t) are bounded. 

Supposing h(t) always positive during the evolution, let us consider separately the 
cases W > and W < 0. In the first case h decreases, since (^]) and (Q) imply 

(14) h = -W-v 2 , 

and so h also is bounded. On the other side, W < implies h 2 < | and so h 
is bounded again. In any case the solution lives in a compact set, and this means 
that sup I = +oo, and it is easily proven that it converges to an equilibrium point 
(0*, 0, h*) of the system. 

Translated into the original formulation, we deduce that if an originally expanding 
solution does not recollapse, then it is regular for all times, the scalar <\> converges to 
a critical point 0* of the potential V, with nonnegative critical value, and — > 0. 

An example is sketched in Figure [j], when the phase portrait for the system @- 
(§) is sketched, for a potential that has two local minima, one with positive critical 
value, the other one with negative critical value. In the first case, every solution 
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starting nearby the critical point (0*, 0, a/2^(0*)/3) approaches it, regardless of the 
curvature (k = — 1, 0, or 1). In the second case, the negative minimum determines 
an intersection between the two branches of the set W = 0, that acts as a "tunnel" 
allowing solution such that W > and starting with h > to recollapse. 

The recollapse problem is an interesting feature of these cosmological models that 
deserves to be further investigated (see Section At the moment, we observe that 
(^|) - but a simple inspection of (Q) also - implies that solutions with k — — 1 may 
recollapse only if they violate the w.e.c, and, by the aforesaid, that occurs only if 
V can attain negative values. We will turn back on this situation later when we will 
examine stability of anti de Sitter model. 

Now we examine stability of the equilibria (0*, 0, h*) as expanding solutions, where 

0* is a critical point of the functional and = \j 2V( ^*^ ■ Recall that, by the assump- 
tions made on V(4>), 0* is either a minimum, or a non degenerate maximum. Let us 
now briefly review the second situation: here, one can simply linearize the system 
(§|)-(|8]) to find that eigenvalues are always real, and one of them is always positive, 
and one always negative, which results in the existence of a 1- or 2-dimensional un- 
stable manifold at the equilibrium point, and then, up to a zero-measure set of initial 
data, the scalar field do not approach the maximum. 

Now, consider the case when 0* is local minimum for V, with non negative critical 
value. If 0* is a non degenerate critical value, asymptotic stability of the equilibrium is 
found simply using the linearized flow of and therefore solutions with initial 

data near the equilibrium point approach it, regardless the case k under consideration. 
Here, we want to extend to the case 0* is possibly degenerate. 

Let us begin considering the situation V ((/)*) > 0, in the open topologies W > 
0. Let us consider first the case V(0*) is strictly positive. We call V > V((p*) a 
regular value such that 0* is the only critical point in the connected component of 
V^O — oo, V]) containing 0*, and let us consider a solution ((j)(t),v(t), h(t)) with 
initial data such that vl + 2V(0 O ) < 2V, and let W = W(0) > 0. We define Q to 
be the connected component of the set {(0, v, h) : W(4>, v, h) E [0, W), e(0, v) E 

[2l / (0 Jt ), 2V}} containing the equilibrium point (0*, 0, \J 2V( f*^ )- It is easy to see, 
using ( |r0| ) and (p"5|), that Vt is compact, and positively invariant by the flow. So, let us 
consider a solution with initial data in Vt. Applying LaSalle's invariance theorem P^ ] 
to the functions W and e, we see that the solution must be such that hW — > and 
hv 2 — > 0. But h is strictly bounded away from zero in fi, then both W and v must 
go to zero, which means that the solution approaches the equilibrium point. Since the 
critical value is strictly positive, then ( |TT1 ) says that the scale factor a{t) diverges at 

infinity like exp(t-y/|F(0*)). 

If V{4> Jf ) = 0, the above argument can be easily adapted. In this case the set 
{(<f),v,h) : W((f),v,h) E [0, W], e(0, v) E [2 V(0*), 2 V}} is connected, and we 
choose f2 to be its subset characterized by the property h > 0. The only point in Q 
with h = is exactly the equilibrium point, and so if h(t) — > (recall h is monotone) 
the solution is forced to approach the equilibrium. If by contradiction h{t) had a 
strictly positive limit, we can argue as before to find v — ► and W — > and so 
h — > 0. Moreover, when the minimum is non degenerate, then the couple (0, v) has 
an oscillatory behavior. This fact was already pointed out by Rendall for the flat case 
[^3p, but it actually takes place also when k — — 1. Indeed, taking 0* = for sake 
of simplicity, and using the variable change = A _1//2 r cos 9, v = rsin#, where 
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Figure 2. Phase portrait sketch when V(4>) = x 2 + (2x 3 ) /3 + x 4 /8. 
The behavior of two curves with initial data in the half space {h > 
0}, living in the regions where W is positive and null, respectively, 
is represented. The two curves both approach the equilibrium point 
(0, 0, 0), where 0* = is the local minimum of V, with zero critical 
value. Both curves spin around the equilibrium point for late times. 

2V(4>) = X(f) 2 + U{4>), with U{4>) = o(0 2 ), then the equation for 6 reads 

= - -/isin20 - —U'a-WrcoaO), 
2 2r v ; ' 

and since both r and h go to zero, we deduce that 9{t) = —y/Xt + o(l/t), and this 
implies the oscillatory behavior of <p (and v). An example of this is sketched in Figure 

I 

Then, it remains proven that every solution of flat or negatively curved Robertson- 
Walker cosmology with initial data near the local minimum of the potential expands 
forever approaching this minimum, and the velocity of the scalar field vanishes at 
infinity. 

If W < 0, we must adapt the above argument, since now we do not have a priori 
information on monotony of h (see (0)). As before, we start from the case V^(0*) > 
0, and take V and initial data such that v 2 + 2V{4>q) < 2V, whereas now W will 
be a negative value to determine in such a way that it acts as a lower bound for W. 
Supposing that initial data are taken near the equilibrium, then h starts positive; since 
W(0) > W, qU^and Q imply that W > W and e(t) < 2V. Last relation, in 
particular, shows that V((p(t) < V and this means that <p{t) is forced to stay in the 
potential well V(<f>*) < V(<p(t)) < V. Then 

h 2 > e - + w>?^A + w = hl + w, 
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that is 

h > I 4/1 + \K =: h. 




Then, chosen sufficiently small such that /i 2 > —W, we have that h{t) > =>■ 
> h > 0. Therefore, we can apply again LaSalle's invariance theorem as in the 
case W > and exploit the fact that h is bounded away from zero. 

All in all, a result of stability of local minima of V(4>) with positive critical value 
holds also in the positively curved Robertson-Walker cosmologies. This result, how- 
ever, cannot be extended as before to the case V(4>*) = 0, because the equilibrium is 
such that = 0, and so nearby solutions may recollapse (i.e. h can change sign). 

4. RECOLLAPSING SOLUTIONS AND THE INSTABILITY OF ANTI DESlTTER 

MODELS 

In this section we give some answer to the recollapsing solutions problem sketched 
before, specifically in the open topology cosmologies, and address the question of 
stability of anti de Sitter solutions with respect to perturbations in the Robertson- 
Walker cosmologies. 

Let us consider a solution of (|6|)-(|8|) with h(0) > 0. Since the sign of W is invariant 
along the flow, we know that W(0) > implies, using (|]), that h{t) is monotone 
decreasing. If the potential V(4>) is nonnegative, we already know from previous 
section that no recollapse will take place, and the solution will expand for all times. 
So, let V be such that V _1 (] — 00, V)) has a connected component containing only a 
critical point <p m which is a local minimum with negative critical value. Taking initial 
data such that v 2 + 2V((po) <2V, then e and W will decrease until h remains positive, 
and the solution is forced to remain in the (compact) set such that W E [0, W(0)], 
e < 2V. But there are no equilibrium points to approach in this set, which means that 
there must be a time T when h(T) = 0, and the solution recollapse (notice that this 
happens violating w.e.c). 

Among these solutions, we can consider anti deSitter models, which occur when 
(p(t) = <f> m , and v(t) = 0. Let us call /i 2 = — |V(0 TO ); then ® takes the form 

h = —h 2 — /i 2 that integrates to give h(t) = —(j, tan(//(£ — c)), where c is the constant 
of integration. Then h(t) diverges to —00 in a finite time t s > 0, and so by <Q) 
W — ► +00, and by dT2|), a(t) goes to zero, but this is actually a "false" singularity, 
that does not correspond to a divergence of stress energy tensor which instead remains 
finite, since v and (j) are constant. 

In order to investigate stability of these solutions, we choose initial data (4>o,vo) 
near (0 m , 0), and suppose that ((f>(t),v(t), h(t)) is a solution of (|6|)-([l3l), defined in an 
interval I, such that <fr(t) — > 4> m , v(t) — ► 0. Then h{t) diverges to —00 as t — » sup I. 
Indeed, by (|8|), and recalling V(4> m ) < 0, h(t) is eventually decreasing, but if h(t) 
admits a finite value it would be bounded, so that sup I = +00 and h(t) would tend 
to a negative value, which cannot happen. This fact implies, again, that a(t) goes to 



zero, that means, using (JTTJ), that 

-i 



lim / h(r) dr = —00. 

t^sup I Jq 



Now let us consider equation (|7]), and call f(t) the (bounded) function —V'((j>(t)). 
Then it must be 

v(t) + 3h(t)v(t) = f(t), 
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which can be integrated to give a diverging solution, once we prove that v ^ 0. 
Indeed, one can consider situations where v = 0, but if v(t) = for all time, (0) 
would give V'((f)(t)) = for all time, that means 4>(t) = <p m , obtaining anti deSitter 
model once again. Then v can be taken nonzero without loss of generality, and this 
results in a diverging solution v(t), which is in contradiction with the assumption 0, v 
bounded. We conclude that anti deSitter solutions are the only solutions of (U)-® 
such that <p goes to a critical point of V(4>) with negative critical value. 

5. Discussion and conclusions 

We analysed here the dynamics of homogeneous Robertson-Walker cosmological 
models with a self-interacting scalar field source within a quite general class of (phys- 
ically valid) potentials V(4>), essentially requiring only the function V to be bounded 
from below and divergent when the field diverges. The analysis has been carried out 
by casting the problem as a regular dynamical system whose solutions describe the 
possible Robertson Walker cosmologies together with the scalar field evolution. 

The asymptotic behavior of the expanding (cosmological) solutions depends cru- 
cially on the sign and the extrema of the function V. First of all, if V is positive, the 
scalar field always approaches the minimum and the solutions do never re-collapse. 
This result extends and completes already existing results, for instance for exponen- 
tial potentials, and holds also for the special case of vanishing V at the minimum if 
the topology is not closed. 

The situation complicates drastically for potentials with a negative lower bound. 
Indeed, solutions which are close to the negative minimum are repulsed; from the 
cosmological point of view, the corresponding universes re-collapse. The unique ex- 
ception is the solutions which "sits" on the negative minimum, which of course is the 
Anti De Sitter space-time. In this sense, AdS turns out to be unstable within homoge- 
neous scalar fields cosmologies. 

The phase space in case of several extrema can be sketched using as an example 
a potential with a positive maximum and one negative and one positive minimum. 
Clearly, from the results above, the phase space contains ever expanding solutions 
approaching the positive minimum, and recollapsing solutions. The two sets are sep- 
arated by solutions for which the scalar field tends to the maximum of the potential. 
It can be easily shown however that the initial data leading to such a situation form a 
two-dimensional surface, and therefore these data, besides being of course unstable, 
are not generic. 
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